Density functional study of oxygen vacancies at the Sn02 surface and subsurface sites 
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Oxygen vacancies at the SnO2(110) and (101) surface and subsurface sites have been studied in 
the framework of density functional theory by using both all-electron Gaussian and pseudopotential 
plane-wave methods. The all-electron calculations have been performed using the B3LYP exchange- 
correlation functional with accurate estimations of energy gaps and density of states. We show 
that bulk oxygen vacancies are responsible for the appearance of a fully occupied flat energy level 
lying at about 1 eV above the top valence band, and an empty level resonant with the conduction 
band. Surface oxygen vacancies strongly modify the surface band structures with the appearance of 
intragap states covering most of the forbidden energy window, or only a small part of it, depending 
on the vacancy depth from the surface. Oxygen vacancies can account for electron affinity variations 
with respect to the stoichiometric surfaces as well. A significant support to the present results is 
found by comparing them to the available experimental data. 

PACS numbers: 61.72.jd, 73.20.At, 73.20.Hb, 71.55.Ht 
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I. INTRODUCTION 

Tin dioxide (Sn02) is one of the most used and inter- 
esting materials for the development of solid state gas 
sensors, transparent conductors, and catalysts.-^ Many 
physical properties of this and other oxides are driven 
by defects which are mostly due to the ease with which 
the oxygen content can be varied. In this respect, Sn02 
surfaces are even more interesting because the presence 
of two possible oxidation states of tin (-1-2 and -1-4), 
which are combined with the reduced atomic coordina- 
tion, favors compositional changes and reconstructions. 
Although the SnO2(110) surface is one of the most stud- 
ied both theoretically^'^''*'^ and experimentally,^'^'^ its 
actual structure is still a matter of debate^ because of 
the controversial aspects concerning defects and oxygen 
adsorption.^ Several reconstructions are possible depend- 
ing on both the preparation conditions and the sample 
history. Nevertheless, in a faceted sample, the lowest in- 
dices (110) and (101) stoichiometric surfaces seem to be 
the most likely, to be favored from the thermodynamic 
point of viewi^ 

Bulk and surface oxygen vacancies are extremely im- 
portant for determining the electrical conductivity of 
tin dioxidei^ Moreover, they are also responsible for a 
very efficient luminescence activity of Sn02 nanobelts. 
First-principles calculations have revealed that bulk oxy- 
gen vacancies and tin interstitials have a low formation 
energyjiS However, few theoretical calculations of the 
energy levels associated to these types of defects have 
been published j-^^^ probably because standard function- 
als used in Density Functional Theory (DFT) fail in giv- 
ing the correct energy gaps .^^'^'^i^^'^^i"'^^ In this paper, we 
report on a first-principles study of the electronic and 
structural properties of Sn02 in the presence of oxygen 
vacancies both in the bulk crystal and at surface and 



subsurface sites, focusing on the effects on the energy 
levels close to the Fermi energy. Accurate calculations of 
the electronic band structure have been done within an 
all-electron approach, which employs the hybrid B3LYP 
exchange-correlation functional. Electron affinities have 
been studied with pseudopotential plane-wave calcula- 
tions. 

In Sec. ini we outline the computational approaches 
used in this work. The structural and electronic proper- 
ties of defected bulk Sn02 are presented in Sec. IIIII The 
role of the oxygen vacancies on the surface properties are 
reported and discussed in the case of the (110) and (101) 
surfaces in Sees. IIVI and FVl respectively. In particular, 
we consider several defect sites which are both at surface 
and in subsurface positions, and analyze their formation 
energy. Electron affinity variations driven by defects are 
studied in Sec. IVII Finally, in Sec. IVIIi we outline some 
conclusions. 



II. COMPUTATIONAL METHODS 

The structural and electronic properties of tin diox- 
ide crystal and surfaces have been studied, from first- 
principles, in the framework of DFT. The first approach 
is based on an all-electron, Gaussian basis set and is im- 
plemented in the CRYSTAL06 (Ref. IT) package. The hy- 
brid B3LYP— li^iiiii^ii^ approximation for the exchange- 
correlation functional has been adopted. This has been 
shown^^'^^ to be extremely valid for band structure cal- 
culations, giving reliable estimations of the energy gaps 
which are in excellent agreement with experimental data. 
We have verified that this approach works well for cas- 
siterite (bulk rutile Sn02), as it will be shown below. 
The all-electron basis set adopted for tin atoms corre- 
sponds to a 3-21G basis setfi^ with the external sp shell 
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FIG. 1: (Color online) Side view of the (a) SnO2(110) and (b) 
SnO2(101) stoichiometric surfaces. The large (small) spheres 
represent tin (oxygen) atoms. The oxygen atoms are labeled 
according to the convention discussed in the text. 



being modified with an exponent a — 0.11. The basis 
set adopted for oxygen atoms is the same as that used 
in Ref. 20, where it was applied to magnesium oxide. 
For the study of the defected structures, the basis set 
was enriched by three sp single Gaussian functions (ex- 
ponents a — 1.0, 0.30, and 0.09) centered at the location 
of the missing oxygen, in order to allow the description 
of trapped electronsi^ 

The second approach, which is based on atomic 
pseudopotentials and a plane-wave basis set as imple- 
mented in the QUANTUM-espresso^^ package, was em- 
ployed due to its being very efficient for treating ex- 
tended, periodic systems. In this case, the calculations 
were performed using the local density approximation 
(LDA) with the Perdew-Wang parametrization for the 
exchange-correlation functional and Bachelet-Hamann- 
Schliiter (for tin) and Vandcrbilt ultrasoft (for oxygen) 
pseudopotentials to represent the ionic cores. We used 
an energy cutoff of 30 Ry for the wave functions and 180 
Ry for the charge density. The need of using a higher 
value of the charge density cutoff is due to the fact that 
we need an accurate estimation of the electrostatic po- 
tentials for the calculation of the electron affinities (see 
below) . 

On the basis of our experience ^^i^ different computa- 
tional schemes converge to comparable results when all 
the main parameters are carefully set. In the present 
case, we checked the charge density cutoff energy and 
the vacuum width parameters, for the plane-wave cal- 
culations, while we carefully set the Gaussian basis set 
by adding ghost orbitals (see below) for the all-electron 
approach. Once verified that both methods give reliable 



results, we used them for different types of calculations. 
The CRYSTAL06 package was used to perform calculations 
of band structures and density of states, thus taking ad- 
vantage of the accuracy of B3LYP exchange-correlation 
functional in the single particle energy calculations. On 
the other hand, the QUANTUM-ESPRESSO package was 
used to perform spatial averages [of potential energy and 
charge density (see below)] over surface planes, align the 
surface electrostatic potential with the bulk one, and, 
therefore, perform calculations of electron affinity varia- 
tions induced by defects. 

The rutile SnO2(110) and (101) stoichiometric sur- 
faces, which are shown in Figs. [Ua) and [Ijb), re- 
spectively, have been modeled with a five layer (seven 
layer) slab in the case of crystalo6 (quantum- 
espresso) calculations. The slab has been studied as 
an isolated structure (two-dimensional periodic struc- 
ture) with the all-electron CRYSTAL06 code, whereas a 
three-dimensional supercell has been considered in us- 
ing QUANTUM-ESPRESSO. In the latter case, we have 
checked that a vacuum width of 6 A between the slab 
replicas is large enough for obtaining converged results. 
For the sampling of the Brillouin zone, we used a 3x3 
Monkhorst mesh (and a 6x6 secondary mesh for evaluat- 
ing the density matrix^^) for the geometrical optimiza- 
tion of the defected structures. Instead, we used a 6x6 
Monkhorst mesh (12x12 secondary mesh) for the band 
structure and density of states calculations. We checked 
that finer meshes do not change the results. 

Defected surfaces have been obtained by removing the 
oxygen atoms in the positions labeled in Fig. [TJ The 
outermost atomic plane of the stoichiometric SnO2(110) 
surface [panel (a)] of Fig. [T] is characterized by rows of 
oxygen in bridging sites between pairs of tin atoms. The 
atomic plane lying just below is made of alternating rows 
of oxygen and tin atoms. The next plane is again made 
of only oxygen, which is again arranged along rowsi^ 
The O atoms belonging to the three atomic planes just 
mentioned will be referred to as bridging, in-plane, and 
sub-bridging, respectively. Bridging, in-plane, and sub- 
bridging O atoms form a primitive structural unit along 
the direction perpendicular to the surface, which will be 
referred to in the following as a surface layer 3^ To label 
these atoms, we use a number for the layer depth (1 is 
the outermost layer) and a symbol for the position within 
the layer: bridging (B), in-plane (IP), and sub-bridging 
(SB). In the same way, the stoichiometric SnO2(101) sur- 
face shown in panel (b) of Fig. [T]is characterized by rows 
of bridging (sub-bridging) oxygen atoms, located above 
(below) planes of tin atoms. The stoichiometric surfaces 
are always favored over other kinds of reconstructions 
when samples are treated in oxygen-rich atmosphere. ^"^ 
In the present paper, we are interested in analyzing how 
oxygen vacancies modify the properties of a stoichiomet- 
ric surface. We calculate the vacancy formation energy as 
E^i = i^vac + ^£^02 - -^st, where E^^^, Eq^, and E^t are, 
respectively, the total energy of the slab containing the 
vacancy, that of an isolated spin polarized O2 molecule. 
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TABLE I: The lattice parameters (a, c), the internal coor- 
dinate u, the energy gap Egap, and the c-axis parallel and 
perpendicular components of the conduction band effective 
masses of cassiterite (bulk rutile Sn02). Theoretical results 
are obtained from the B3LYP, all-electron calculation. Ex- 
perimental data are reported for comparison. The lattice pa- 
rameters a and c are in A, the energy gap is in eV, and the 
masses are in free-electron-mass units. 
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and that of the stoichiometric slab. The present expres- 
sion is equivalent to Eq. (2) of Ref. [l^, when the partic- 
ular case of a neutral oxygen vacancy is considered. 

The geometries of all the systems considered in this 
work have been optimized by fully relaxing the atomic 
positions until the residual forces were less than 10~^ Ry 
/ Bohr. We conclude this section by pointing out that 
the inclusion of the spin polarization does not modify our 
results. Indeed, we checked that, at variance with tita- 
nium dioxide (where in the presence of oxygen vacancies, 
a triplet ground state may arise), tin dioxide surfaces 
always show a spin singlet ground state which is fully 
justifying the use of spin unpolarized calculations. 



III. BULK PHASE 

In Table [H the calculated lattice parameters, energy 
gap, and conduction band effective masses of bulk ru- 
tile Sn02 are compared to the experimental data, and 
show good agreement. The theoretical results, which are 
obtained with the all-electron B3LYP calculations, show 
deviations from the experimental data of only about 1.5% 
for the structural properties and less than 3% for the en- 
ergy gap. 

In Fig. [2l we show the band structure of bulk Sn02. 
As already known, Sn02 is a direct band gap semicon- 
ductor, with the lowest energy electron-hole transition 
occurring at the F point. It is interesting to compare the 
present results with previous band structure calculations. 
Former tight binding calculations^^ fairly reproduced the 
lowest bands (located at 15-20 eV under the top valence 
band) and the overall shape of the conduction band. Yet, 
they fail in the calculation of the top valence band, due 
to a wrong estimation of the charge transfer toward the 
oxygen atoms. This failure is critical in the presence 
of an oxygen vacancy, mostly if any geometrical recon- 
struction at the vacancy site is neglected.— >^ Other inter- 
esting results emerge from a comparison with the GGA 
density-functional calculations of Ref. ISOl . In that paper, 
Maki-Jaskari and Rantalcti^ showed that ultrasoft pseu- 




FIG. 2: (Color online) The band structure of bulk Sn02 ob- 
tained from the all-electron, B3LYP calculation. Blue solid 
(red dashed) lines refer to occupied (unoccupied) states. A 
dotted line indicates the highest occupied level. 



dopotentials (USP) provide more accurate results than 
their norm- conserving (PSP) counterparts, in the energy 
range close to the top valence band. We confirm their 
conclusion since the shape of the present B3LYP band 
structure (see Fig. [2]) is in good agreement with the USP 
calculation within the full valence band spectrum. 

In Fig. [21 we show the calculated and measured den- 
sity of states (DOS) of bulk Sn02. The theoretical re- 
sult [curve (a)] has also been broadened using a Gaus- 
sian distribution with a half width at half maximum of 
0.3 eV [curve (b)] to facilitate the comparison with the 
experimental^ data [curve (c)]. Figure [3] shows good 
agreement between theory and experiment. In partic- 
ular, we find a valence band width of 8.1 eV in good 
agreement with both experimental data (7.5 eV reported 
in Ref. 0) and previous first-principles calculations (7.9 
eV and 8.8 eV reported in Ref. M using PSP and USP, 
respectively) . 

Before discussing the results concerning the surface 
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FIG. 3: DOS of bulk tin dioxide, (a) Present B3LYP cal- 
culation. The top of the valence band has been chosen as a 
reference energy, (b) The DOS is convoluted with a Gaus- 
sian curve whose half width at half maximum is 0.3 eV. (c) 
Experimental data taken from Ref. S- 
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FIG. 4: Band structure (left panel) and DOS (right panel) 
of bulk Sn02 with an oxygen vacancy in a bridging site. A 
2x2x3 supercell has been used in the calculation. The shaded 
regions represent the projected band structure and DOS of 
the undefected structure. 



band structures with and without oxygen vacancies, it 
is worth addressing their role in bulk Sn02. Isolated va- 
cancies have been modeled using a 2x2x3 (almost cubic) 
supercell, which contains 12 primitive unit cells. From 
each supercell, we remove a bridging oxygen atom, so 
that the final structure results in a vacancy-vacancy dis- 
tance of about 9.5 A. We computed a vacancy forma- 
tion energy of 3.42 eV for this configuration. Figure [4] 
shows the band structure (left panel) and the DOS (right 
panel) in the near-gap region for the defected structure 
(lines). The shaded regions represent the projected band 
structure and the DOS of the ideal crystal, respectively. 
All the energy levels located at negative (positive) ener- 
gies are occupied (unoccupied) at a temperature of °K. 
The uppermost occupied state of the defected structure 
is taken as reference energy. The comparison between the 
pure and the defected crystals leads to important conclu- 
sions. First, vacancies are responsible for the appearance 
of an almost flat occupied band which is located about 
1 eV above the top valence band of the ideal crystal. 
The corresponding peak in the DOS is well visible in the 
right panel of Fig. [U which shows a narrow well-defined 
distribution. Interesting results also emerge for the unoc- 
cupied states. Indeed, Fig. [4] evidences the appearance 
of states related to the oxygen vacancy and resonating 
with the conduction band. The coupling between the 
conduction band and the resonance level localized at the 
vacancy site pushes the bottom of the conduction band 
at a higher energy. A clear evidence of the effects of the 
vacancy on the energy levels emerges from the DOS. The 
vacancy is responsible for both the blue-shift of the den- 
sity of unoccupied states threshold and the appearance 
of a strong peak at about 1.1 eV above this threshold. 



IV. SNO2(110) SURFACE 

As outlined in sec. [Ill different vacancy sites (bridg- 
ing, in-plane, and sub-bridging) are possible for the 
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FIG. 5: (Color online) Band structure of the SnO2(110) stoi- 
chiometric surface. The reference energy is the vacuum level. 
The last occupied state is at -8.4 eV and is indicated with 
a dotted line. Occupied (unoccupied) bands are reported as 
blue solid (red dashed) lines. The shaded region corresponds 
to the projected bulk band structure. The alignment of the 
projected band structure is fixed to match the Is core levels of 
the Sn atoms in bulklike positions. The bands are calculated 
using the B3LYP scheme. 



SnO2(110) surface. We have considered all these vacancy 
positions, at several depths inside the slab. In all cases, a 
2x1 surface reconstruction (the x and y axes have been 
chosen along the [001] and [110] crystallographic axes, 
respectively) is considered. 

The band structure of the stoichiometric (undefected) 
tin dioxide surface is reported in Fig. [5l The reference 
energy is the vacuum level. The comparison with the 
bulk Sn02 band structure evidences that surface bands 
show up inside the forbidden band gap. The projected 
DOS (PDOS) onto atomic orbitals shows that the top 
valence band mostly consists of p states localized on the 
external bridging (IB) oxygen atoms. The bottom of the 
conduction band mainly consists of s orbitals localized 
on the in-plane tin atoms. The formation of the intragap 
surface states leads to a surface band gap much smaller 
than the bulk one. We find that the stoichiometric sur- 
face has a band gap at the F point of 2.41 eV. Figure 
[5] is in qualitative agreement with previous theoretical 
calculations,'^'^ with the difference that the present ap- 
proach gives reliable results for the conduction band en- 
ergies. An indication of the good quality of the present 
calculation can be found in the position of the lowest 
conduction band with respect to the vacuum level. We 
find a value of -5.7 eV [-5.5 eV for the SnO2(101) sur- 
face], whereas, secondary electron cutoff measurements 
in photoemission gave a work function value (defined as 
the difference between the vacuum and the Fermi level 
energies) of 5.7±0.2 eV for the stoichiometric SnO2(101) 



o 
Pi 



w 







FIG. 6: (Color online) Band structure of the SnO2(110) sur- 
face with an oxygen vacancy at the outermost bridging (IB) 
(left panel) or in-plane (IIP) (right panel) positions. The ref- 
erence energy is the vacuum level. The highest occupied states 
he at -5.74 eV (left panel) and -7.35 eV (right panel) and are 
indicated by dotted lines. Occupied (unoccupied) bands are 
reported in blue solid (red dashed) color. The shaded region 
corresponds to the projected bulk band structure. The align- 
ment of the projected band structure is fixed to match the Is 
core levels of the Sn atoms in bulklike positions. The bands 
are calculated using the B3LYP scheme. 



surface.-^ According to the literature, the Fermi level lies 
just below the lowest conduction level, ^ (tin dioxide be- 
haves as a n-type material), and thus, the agreement with 
the present calculation is good. 

In Fig. [21 we show the band structures of a de- 
fected SnO2(110) surface, where either a bridging (IB, 
left panel) or an in-plane (IIP, right panel) oxygen atom 
has been removed. There are substantial differences with 
respect to the stoichiometric surface, especially concern- 
ing the intragap states. As a result of the removal of the 
most exposed (IB) bridging oxygen atoms, a dispersed 
band arises within the forbidden gap. The trend found 
here is in agreement with the experimental literature, 
reporting that, upon reducing a stoichiometric surface, 
energy levels extend from the top bulk valence band up 
to about 0.5 eV below the conduction band minimum.^ 
We found a different picture when the vacancy lies in the 
subsurface positions. In this case, the levels within the 
gap more closely resemble those of bulk Sn02, with a lit- 
tle dispersed band placed at about 1 eV above the top 
bulk valence band, as shown in the right panel of Fig. [S] 

In Fig. [71 we report the vacancy formation energy 
[panel (a)], the energy gap [panel (b)], and the highest 
occupied state (HOS) binding energy calculated with re- 
spect to the vacuum level [panel (c)] , upon changing the 
depth of the oxygen vacancy inside the slab. Squares, cir- 
cles, and triangles correspond to the bridging, in-plane 
and sub-bridging oxygen vacancies, respectively. The 
black full symbols show calculations performed using the 
all-electron Gaussian based CRYSTAL06 package, with the 



B3LYP exchange correlation functional. As a compari- 
son, the red empty symbols represent calculations per- 
formed using the plane- wave pseudopotential quantum- 
ESPRESSO package, with a LDA for the exchange correla- 
tion functional. Interesting results can be inferred from 
this figure. First, the vacancy formation energy [Fig. 
[Tlja)] changes with the depth going toward the bulk limit 
of 3.42 eV for sufficiently deep vacancy sites. Although 
there is some scattering in the data, there is always a 
monotonic increasing behavior of the curve at a fixed va- 
cancy position [compare 1B-2B, 1IP-2IP, and 1SB-2SB 
in Fig. [7fa)]. In other words, for a given position, the 
vacancy formation energy either increases or remains con- 
stant. This result is at variance with that found in the 
case of titanium dioxide, where a more oscillating behav- 
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The results on tin dioxide vacancy formation energies 
already available in the literature span a wide range of 
values. There is a strong dependence on the vacancy 
type, depth, and concentration, but also on the com- 
putational method used for the calculations. For the 
outermost bridging oxygen vacancy (IB in Fig. [1}, the 
plane- wave calculations of Oviedo and Gillan^ and Maki- 
Jaskari and Rantala^ gave values of 2.6 eV and 3.0 eV, 
respectively. On the other hand, the localized orbital 
calculations of Sensato et alM- lead to a value of 4 eV 
for the formation energy of a fully reduced surface, in 
which all bridging atoms had been removed. A possible 
source of discrepancy among the methods can be found 
in the erroneous evaluation, within a Gaussian basis set, 
of the charge density in the space region of the missing 
oxygeuj^ due to the non negligible ionicity of tin dioxide. 
In order to fix this error, we added three Gaussian type 
functions to the original basis set centered at the vacancy 
site (ghost orbitals)j^ We verified that the addition of 
ghost orbitals reduces the vacancy formation energy (it 
is 0.5 eV for the IB vacancy), which leads to a better 
agreement with the plane-wave results. Another impor- 
tant source of discrepancy lies in the use of either PSPs 
or USPs in the plane-wave approach, which causes strong 
variations of the vacancy formation energiesj^ The data 
reported in Fig. [3 show the overall fair agreement that 
we obtained for the vacancy formation energy using ei- 
ther the all-electron B3LYP Gaussian basis set or the 
pseudopotential plane-wave LDA calculations. 

Figure ^h) shows the energy gap i^gap as a function 
of the vacancy depth. i?gap strongly depends on the hy- 
bridization of atomic orbitals surrounding the vacancy, 
which explains the dependence on the vacancy site. An 
interesting point is that, with the exception of the more 
exposed (IB) oxygen whose removal leads to the lowest 
value of the energy gap (due to the formation of many 
states within the gap), for each vacancy type there is 
a negligible dependence on the depth. There is a gap 
for bridging, in-plane, and sub-bridging vacancies. Since 
such a behavior seems very promising for technological 
applications, further calculations are being performed on 
deeper vacancy sites. Figure [3[c) shows the HOS bind- 
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FIG. 7: (Color online) (a) Vacancy formation energy, (b) en- 
ergy gap, and (c) highest occupied state position with respect 
to the vacuum level. Bridging (squares), in-plane (circles), 
and sub-bridging (triangles) oxygen vacancies are consid- 
ered. Black full and red empty symbols represent calculations 
performed using either the all-electron B3LYP CRYSTAL06 
package or the pseudopotential plane-wave LDA QUANTUM- 
ESPRESSO package. The lines are guides for the eyes. 



ing energy with respect to the vacuum level. There is an 
overall reduction of the HOS energy upon increasing the 
depth at which the vacancy is formed. This reduction of 
the HOS is in line with the electron affinity increase (see 
below). Unfortunately, in the literature, there is a lack 
of experimental measurements of the HOS energy versus 
the vacancy depth. 

In Fig. [51 we show the DOS of both the stoichiomet- 
ric surface [panel (a)] and of the surface with a bridging 
(IB) oxygen vacancy [panel (b)]. The reference energy 
is the vacuum level. As already mentioned, the surface 
reduction (e.g., after sputtering) leads to the formation 
of intragap surface states, ^-^ but the oxidized (stoichio- 
metric) surface shows surface states as well (see Fig. [SJ. 
This explains the experimental finding that the anneal- 
ing process in air does not fully remove the intragap 
states. We find that the stoichiometric surface is fea- 
tured by occupied intragap states at about 1 eV above 
the top bulk valence band. This is in agreement with 
previous calculation a^^i'^° and experiments.^ Upon reduc- 
tion of the surface, a continuum of occupied states which 
covers a wide energy range within the bulk band gap is 
experimentally reported.^ This clearly emerges from Fig. 
ISl^b). The main peak above the bulk valence band is 
lower in value than the corresponding peak of the sto- 
ichiometric surface. At the same time, a wider energy 
range within the bulk gap is spanned by the occupied 
states. In the figure, the highest occupied and lowest 
unoccupied states are indicated with arrows and the for- 
bidden energy window between those states is reported 
using a dashed line. The reduction in the energy gap 
due to the presence of the vacancy is clearly seen. From 
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FIG. 8: (Color online) DOS of (a) the stoichiometric and 
(b) the defected (IB O vacancy) SnO2(110) surface. The 
(smoother) thick lines are a convolution of the DOS with a 
Gaussian broadening of 0.4 eV. We indicate with the shaded 
region the original (nonconvoluted) DOS. The highest occu- 
pied and lowest unoccupied states are reported in the figure 
using vertical arrows, whereas the forbidden energy window 
is indicated by a horizontal dashed line. 



Fig. [H we also observe the shift to higher energies of 
the whole bulk spectrum of the defected with respect to 
the stoichiometric surface. This is in agreement with the 
decrease in work function by about 1 eV reported in the 
experimental literature i^i*^ 



V. SNO2(101) SURFACE 

For the study of oxygen vacancies at the SnO2(101) 
surface, we have considered a v2 x v2 reconstruction 
(double cell, i.e., the area is doubled with respect to the 
unit cell of the bulk geometry). We studied the vacancy 
formed by removing a top bridging oxygen atom from 
each unit cell with a vacancy-vacancy distance of 7.4 A. 

In Fig. [9l the band structures of the stoichiometric 
(left panel) and the defected (right panel) surfaces are 
shown. In order to compare both band structures, a 
\f2 X a/2 reconstruction is considered for the stoichiomet- 
ric surface as well. Similarly to what is discussed in the 
case of the SnO2(110) surface, the presence of the surface 
leads to the formation of surface bands within the bulk 
band gap. The stoichiometric surface has a direct gap 
at r, with an energy gap of 2.7 eV. The presence of the 
oxygen vacancy leads to the formation of a doubly degen- 
erate occupied energy level, which is well separated from 
all the other bands. This level shows a parabolic dis- 
persion near F, but it is almost flat elsewhere inside the 
first Brillouin zone.. From the PDOS analysis, we found 
that the corresponding states are well localized around 
the vacancy site, which explains in this way the lack of 
dispersion. The oxygen vacancy is also responsible for 
a level resonant with the conduction band. At variance 
with the stoichiometric surface, the bottom of the con- 
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FIG. 9: (Color online) Band structure of the stoichio- 
metric (left) and defected (bridging oxygen vacancy, right) 
SnO2(101) surfaces. The last occupied states are at -8.2 eV 
and -6.76 eV, respectively, and are indicated by a dotted line. 
Occupied (unoccupied) bands are reported in blue solid (red 
dashed) color. The shaded region corresponds to the pro- 
jected bulk band structure. The alignment of the projected 
band structure is fixed to match the Is core levels of the Sn 
atoms in bulklike positions. The bands are calculated using 
the B3LYP scheme. 



duction band for the defected surface lies at the X point. 
The first allowed transition for the SnO2(101) surface 
when a top bridging oxygen vacancy is formed is at the 
X point. Interestingly, the electron and hole states in- 
volved in this direct-gap recombination are quite local- 
ized around the vacancy site. This localization leads to 
a strong oscillator strength and may explain the photo- 
luminescence activity of tin dioxide nanowiresi^^^ 



VI. ELECTRON AFFINITY 

Electron affinity (EA) variations induced by the pres- 
ence of vacancies have been studied for oxygen vacancies 
located on the first two layers of the SnO2(110) surface 
[see Fig. [IJa)]. This part of the work has been done us- 
ing the QUANTUM-ESPRESSO code. The EA can be easily 
calculated once the vacuum level has been determined 
from the self-consistent electrostatic potential)^ The re- 
sults are collected in Table HIl In this table, we give, for 
each vacancy type, the depth of the removed oxygen cal- 
culated from the top bridging oxygen. It is of interest to 
see how the different vacancies contribute to the EA. 

The top panel of Fig. [Tni shows the planar average 
of the electrostatic potential. The averages are taken 
on planes parallel to the surface (the z coordinate lies 
along the surface normal and z~Q corresponds to the slab 
center). The potential is shown for the stoichiometric 
case and for the ISB and 2B vacancies. If we take the 
stoichiometric surface as a reference, it is seen that the 
2B vacancy induces an EA variation of 0.62 eV, whereas 
the variation due to the ISB vacancy is of -0.57 eV. This 



TABLE II: EA and its variation with respect to the stoichio- 
metric surface (AEA), calculated for the SnO2(110) surface 
with an oxygen vacancy, whose depth is reported in the sec- 
ond column. The vacancies are labeled according to Fig. [IJa) . 
All the energies are measured in eV, and the depth is in A. 



Vacancy type 


Depth 


EA 


AEA 


IB 


0.0 


3.69 


-0.90 


IIP 


1.23 


4.00 


-0.59 


ISB 


2.70 


4.02 


-0.57 


2B 


3.60 


5.20 


0.62 


2IP 


4.87 


4.54 


-0.05 


2SB 


6.30 


4.91 


0.33 


Stoichiometric 




4.59 


0.0 



difference in both the amplitude and sign arises from the 
different electron charge density and atomic relaxations 
following the oxygen removal. To better illustrate this 
point, we plot in the bottom panel of Fig. [10] the planar 
averaged electron charge density. An electron density 
depletion occurs in correspondence of the vacancy site, as 
indicated by the vertical arrows in the figure. Moreover, 
the atomic relaxation can be inferred from the shift of 
the charge density peaks. The combination of both the 
charge depletion and the atomic relaxation gives rise to 
an electric dipole variation (that we have calculated as in 
Ref. [li) AP = 0.27 a.u. for ISB and AP = -0.29 a.u. 
for 2B. The change in the electron affinity due to a change 
AP in the surface dipole is Ax — —AttAP/A, where A 
is the surface unit cell area (all quantities measured in 
atomic units). '^^'^^ Within the surface dipole model, we 
therefore obtain Ax = —0.57 eV for the ISB vacancy 
and Ax = 0.61 eV for the 2B vacancy with remarkable 
agreement with the data of Table [TTl 

Few experimental works have been published on the 
dependence of electron affinity and work function on the 
surface oxidationi^i^ In particular, in Ref. ^, a care- 
ful study was done on the effects of the annealing tem- 
perature on the electron affinity and conductivity of an 
oxidized (presumably stoichiometric) SnO2(110) surface. 
Cox et al^ found that upon increasing the annealing 
temperature, the amount of oxygen lost by the surface 
increases. They concluded that for temperatures up to 
about 700 °K, only the bridging oxygen (IB in Fig. 1) 
is removed from the surface followed, at higher tempera- 
tures, by the in-plane oxygen. This trend is in qualitative 
agreement with the vacancy formation energies shown in 
Fig. [TJa). The results of Ref. ^ show that the varia- 
tion of electron affinity with respect to the stoichiometric 
surface is negative, with a maximum variation of about 
0.4 — 0.5 eV upon increasing the annealing temperature. 
This is compatible with our results in Table ITIl However, 
it should be stressed that a strict comparison between 
our results and the experimental data is made difficult 
by both the uncertainties in the actual surface structure 
and the computational difficulties in analyzing a large 
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FIG. 10: (Color online) Planar average of the self-consistent 
electrostatic potential energy (upper panel) and electron den- 
sity (lower panel) for the stoichiometric SnO2(110) surface. 
The ideal surface is compared with the defected one (ISB 
and 2B oxygen vacancies). All the potentials calculated for 
the studied configurations are aligned with the bulk potential 
in the middle layer of the slab. The energy reference is taken 
on the top of the bulk Sn02 valence state. The electron affin- 
ity X is defined as the difference between the vacuum level 
and the bulk conduction state energy. The position of the 
oxygen vacancy within the slab is indicated with an arrow 
in the lower panel. The calculations have been done using a 
pseudopotential plane- wave approach. 



number of different vacancy configurations. 

Finally, we conclude by noticing that the recently re- 
ported decrease in the work function of about 1 eV upon 
reducing the surface^liH is in agreement with the results 
shown in Table |IT1 where the electron affinity decreases 
by 0.9 eV when an outermost oxygen atom is removed 
from the stoichiometric surface. 



SnO2(110) and (101) surfaces and their role on the sur- 
face electronic properties. As expected, the formation 
energy tends to increase, when the vacancy moves from 
the surface to subsurface layers. The energy gap, instead, 
is very sensitive to the vacancy geometry (and reconstruc- 
tion), but it only weakly depends on the layer depth in- 
side the material. Interesting results on the alignment of 
energy levels and the appearance of intragap states have 
been discussed in connection with the presence of oxygen 
vacancies. 

Some results have also been obtained concerning the 
role of oxygen vacancies in bulk tin dioxide. We have 
found that oxygen vacancies do not lead to the forma- 
tion of occupied donor levels which are just below the 
conduction band. Instead, the present calculations have 
shown that occupied levels associated to the vacancy ap- 
pear at about 1 eV above the top valence band. 

Another important conclusion of the present work is 
that the process of surface reduction, which leads to the 
removal of the most external oxygen atoms, pushes the 
density of states spectra to higher energies, as already 
known from experimental results. 

Finally, our calculations have shown that oxygen va- 
cancies may significantly modify the material electron 
affinity. Interestingly, subsurface vacancies can either in- 
crease or decrease the electron affinity and this effect can 
be described in terms of a surface dipole layer that takes 
contributions from both the atomic relaxation and elec- 
tron charge depletion around the missing oxygen. These 
findings could help in understanding the basic physics of 
the Sn02 surfaces. 
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